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Abstract 

The energetic cost of creating a defect within a host material is given by the formation energy. 
Here we present a formulation allowing the calculation of formation energies in one-dimensional 
nanostructures, which overcomes the difficulties involved in defining the chemical potential of the 
constituent species and the possible passivation of the surface. We also develop a formula for 
the Madelung correction for general dielectric tensors and computational cell shapes. We apply 
this formalism to the formation energies of charged Al impurities in silicon nanowires, obtaining 
concentrations significantly larger than in their bulk counterparts. 

PACS numbers: 68.35.Dv,73.20.Hb,61.46.Km 
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The study of the energetics of the formation of defects is a very well-developed topic in 



bulk semiconductor physics 



Q,y,Q. 



Ill ], but it is considerably less mature in nanostructures. 



A sensible definition of the chemical potential, the possible presence of passivating agents 
on the surface and a proper treatment of the defect charge state are some of the issues that 
prevent the extension of the standard theory from being straightforward. The formation 
energy of a defect Q] is a quantity of paramount technological importance; it determines the 
structural configuration and the charge state that a given set of impurities will favor, and it 



is used in the computation of impurity equilibrium concentrations 



5], solubilities 



0,3, 



diffusivities J3, |sj], dopant compensation mechanisms jsj], etc. In addition, the calculation 
of the formation energy is required whenever a comparison between configurations with 
different number of atoms/chemical species is wanted. 

In bulk host materials, formation energies are calculated according to the well-established 
theory due to Zhang and Northrup ]3], where they are formulated in terms of the chemical 
potentials of the constituent species and the total energy of the system with the impurities. 
On the other hand, for one-dimensional (ID) semiconductor systems, the nonequivalence of 
the different constituent atoms in, say, a silicon nanowire (SiNW), in addition to surface 
passivation issues, render the straightforward application of the Zhang-Northrup formalism 
troublesome. In particular, the definition of the chemical potential of the atomic species in- 
volved is ill-defined. This difficulty is considerably lessened in the study of nanotubes, where 



lis has allowed initial calculations of defect forma- 



each atom is equivalent to the others. T 
tion energies in C and BN nanotubes [12 



131 ]. but analogous calculations for semiconductor 
nanowires have been missing so far. Additionally, the most stable configuration of a defect 
in a semiconductor may have a charge state different from zero, depending on the doping 
condition of the material. In a periodic boundary condition (PBC) formalism, a finite net 
charge in the simulation cell would give rise to a divergent Coulomb energy because of the 
interaction with its periodic images. While the correct procedure for the removal of this 
contribution to the total energy is still a matter of debate, recent reports indicate that the 
uniform background charge 14. 15] and the local-moment counter charge yield similar 



results 



171 ] for bulk materials. However, a treatment for charged defects in one-dimensional 



systems is lacking so far. 

In this Letter we propose a framework for the calculation of formation energies of neutral 
and charged point defects in ID systems. As a case study we discuss the formation of Al point 
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defects, which can^ provide ^-doping and can be found as contaminant from Al-catalyzed 
growth process 



231 ]. thus conveying a considerable technological interest. Specifically, we 
will deal with substitutional and interstitial defects at different radial positions in (110) and 
(111) SiNWs of 1.0 and 1.5 nm diameter, identifying whether there is a tendency to surface 
segregation and their most stable charge state for different doping conditions. 

Formation energy in bulk - According to the well-established theory due to Zhang and 
Northrup [3J, the formation energy of a charged defect in a semiconductor AE? is given by 

AEf = E° - E\T - £ ntfn + q(e v + // e ), (1) 

where Ef ot (Ef^ t re ) is the total energy of the defective (clean) system, ni is positive (negative) 
for atoms added to (removed from) the clean system, //, is the chemical potential of the 
reservoir supplying the impurities, e v is the top of the valence band of the clean host and 
He is the chemical potential for electrons. 

Chemical potential - Extending the use of this formalism to one-dimensional nanostruc- 
tures involves some subtleties related to the definitions of the chemical potential, especially 
in the case of vacancies, substitutionals and self-interstitials. When a vacancy is formed in 
bulk, the removed atom is implicitly assumed to be added to the crystal; thus the use of the 
bulk chemical potential is justified. For a NW, the displaced atom can be added to multiple 
nonequivalent positions, as opposed to bulk, where all lattice sites are equivalent. This in- 
duces an ambiguity with respect to the choice of the chemical potential, but even if somehow 
a preferred site for the displaced atom were selected, the calculation of its contribution to 
the total energy would be ill-defined from a computational point of view. In addition, if the 
NW is passivated, it is not clear how many passivating atoms or fraction thereof should be 
assigned to the displaced atom. Clearly, any attempt to use Eq. (jTJ) directly will not be able 
to deal with the possibility that the NWs may be surface passivated. 

We can circumvent these issues by considering the formation of n^f defects instead of 
a single defect (see Fig. [1]), where n^f is the number of silicon atoms in a NW primitive 
cell — for the sake of simplicity we will refer to Si atoms, meaning any host material. These 
displaced atoms are then added at the end of the NW to form an extra primitive cell with the 
needed passivating agents. We notice that this is equivalent to defining a chemical potential 
of the wire primitive cell. 

With the above or analogous constructions, it is easy to see that a general expression for 
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(a) (b) 




FIG. 1: Construction for dealing with defects in ID semiconductor structures, (a) Vacancies and 
substitutionals: we create as many defects as necessary for adding an extra primitive cell to the 
wire, (b) Self-inter stitials: we create as many defects as necessary for removing a primitive cell 
from the wire. 

the formation energy of a given type of defect will be 
AE' = Eg t - NE P N C W - Y, nxf^x- 

^(E^-n p H C fi H ) + q(e v + [x e ), (2) 

n Si 

where N is the number of primitive cells used in the clean system calculation, E^yy is the 
energy of the NW clean primitive cell, satisfying Ef^" e = NE^w, n x is the number of non- 
host atoms of species X added to the clean system, /i x is their chemical potential, risi is 
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the number of Si atoms involved in the defect formation (e.g. +1 for a self- interstitial), n^f 
is the number of passivating atoms in a NW primitive cell (taken to be hydrogens) and hh 
is the corresponding chemical potential. All quantities appearing in Eq. (j2J) are now well 
defined and easily extracted from a total energy calculation. Note as well that, in the bulk 
limit, rigf — > oo and Eq. ([T]) is recovered. 

Charged defects - The study of charged defects in one-dimensional nanostructures in 
a PBC formalism must overcome various particularities not present in the bulk case, all 
arising from the fact that a dielectric tensor e will be needed for the correct description of 
the interaction between the different instances of the charged defect. 

The usual procedure for dealing with these effects in bulk materials consists in using a 
neutralizing jellium background to recover the charge neutrality condition, and then correct 
a posteriori for the spurious terms arising in the total energy by means of a Madelung 



correction [18j divided by the value of the (isotropic) macroscopic dielectric constant of the 
host material [jjj]. 

In the case of a NW, the numeric value of the Madelung constant will depend on the 
relation between the lattice parameters and the chosen e, and thus it cannot be looked 
up in tables. Starting from the solution to the Poisson equation within a homogeneous, 



anisotropic medium [19|, and following analogously to the standard procedure 18], one can 
easily obtain the following expression for the Madelung constant in the general case 



El erfc(7A/Ri ■ e _1 • Rj 



R Vdete VRi ■ e' 1 ■ R; 

\- 4tt exp(G^ ■ 6 ■ G^/47 2 ) 2 7 tt 

V V c Gi e Gi vVdei7 K7 2 ' 1 ' 

where the sum over R; (Gj) extends over all vectors of the direct (reciprocal) lattice except 
for zero, 7 is a suitably chosen convergence factor and V c is the volume of the primitive cell. 
This approach, besides dealing with generic dielectric tensors, allows us also to easily tackle 
non- conventional cell shapes. This is a very common situation in ID systems, where the 
axial lattice parameter obeys the periodicity of the crystal structure, while the transverse 
dimensions are normally much larger. This allows a proper buffer vacuum to avoid spurious 
interactions with the image neighboring system [^(J . Concerning the choice of the dielectric 
tensor e, one would be tempted to use, for example, the modified Penn model 



2 if ] in order to 



obtain an approximation to the value of the dielectric constant for directions perpendicular to 



FIG. 2: Cross section view of the 1.0 nm (110) (left) and (111) SiNWs (right) studied. The 
substitutional and interstitial Al defects considered are shown in light orange spheres (blue and 
white spheres representing Si and H atoms, respectively). Notice that that the tetragonal interstitial 
is not visible in the (111) SiNWs because it lies along the wire axis and it is covered by host atoms 
and by the innermost substitutional Al defect. 

the growth axis. This would be indeed the correct approach if we were to study the physics 
of excitons in the NW. However, it must be kept in mind that in the case at hand two 
images corresponding to different instances of the NW will interact mainly through vacuum 
if enough buffer space is left. Thus, we will take the dielectric tensor to be diag (1, 1, e r ), 
where e r is the bulk dielectric constant of the constituent material of the NW. 

Al defects in (110) and (111) SiNWs - We apply the formalism described above to the 
study of Al point defects in 1 nm and 1.5 nm SiNWs grown along the (110) and (111) axes. 
Two reasons make Al impurities a very interesting case study: (i) group-Ill elements can 
be efficient p-type dopants for silicon, and the use of Al for doping in nanowires has indeed 



been proposed 



221 ]; (ii) Al has pro ven to be a feasible alternative to Au as a catalyst for 



the epitaxial growth of SiNWs [23|, having the considerable advantage of not introducing 
undesired mid-gap states that can act as traps, and requiring lower growth temperatures. 
We calculate the total energy within density-functional theory, as implemented in the Siesta 
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package [24j], using norm- conserving pseudopotentials, an optimized double-^ basis set [25] 
and the spin-polarized version of the generalized gradient approximation (GGA) 26j] for 
the exchange-correlation energy. The supercell size was chosen to guarantee a separation of 
~23.7 A (~28.3 A) between the impurity and its periodic image along the (110) ((111)) axis, 
while the transverse dimensions were held at 50 A. The atomic positions were relaxed until 
the forces on all the atoms were lower than 0.02 eV/A and the axial lattice parameter was 

n 

optimized for the pristine wire for each growth orientation [27[ . If an Al nanoparticle is used 
to catalyze the growth reaction, the nanoparticle will be the main source of Al contaminants 
in the SiNW. Therefore, we define the Al chemical potential with respect to an Al particle 
of approximately the same diameter of the wire. 

At first we have considered Al substitutional impurities, Alg;, at two different radial 
coordinates, and the tetragonal interstitial, Alf , in 1 nm (110) and (111) SiNWs (see Fig. [2]). 
As expected, Al is a single acceptor, with substitution at a host lattice site being largely 
favored over the interstitial position (see Fig. [3]). Despite the thin diameter of the wire, the 
formation energy of the Als; is hardly sensitive to the change in the radial coordinate of the 
lattice site; the only noticeable change is a slight lowering of the 0/— occupation level of the 
(110) wire when the Al substitutes moves outwards (dashed line in Fig. [3]). We have not 
observed any marked tendency to surface segregation, with the only exception of n-doped 
(110) wires, where the center Als; releases 0.15 eV when moving to more external position. It 
should be noted, however, that Fernandez-Serra et al. showed that the segregation energies 
rapidly increase in presence of surface defects [281 ]. a situation not considered here. The 
most important fact, however, is that the formation energy of the Als; is negative for n- 
type doping condition, meaning that Al is more stable as a point defect in the wire rather 
than in the catalyst nanoparticle. Hence, n-type doping in presence of an Al catalyst might 
be impractical, as it results in easy Al incorporation in the wire, leading to uncontrolled 
compensation. A role is certainly played by the reduced value of the Al chemical potential 
fi ai in the nanoparticle with respect to bulk ct-Al (see Eq. [2]). From the behavior of the 
formation energy vs. the Fermi level, we can infer that choosing properly the substrate 
doping allows the control of the wire dopant concentration. Straightforward application 
of the Zhang-Northrup formalism, using bulk derived chemical potentials, results in an 
underestimation of the formation energies of ~ 0.1 eV. 

In an attempt to address the formation energy vs. the wire diameter, we have studied 
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FIG. 3: Formation energy of Al substitutional (Algi) and Al tetragonal interstitial (AlJ) in (110) 
and (111) 1.0 nm SiNWs. Substitutional defects have been studied at two different radial coordi- 
nates: at the center of the wire (continuous line) and at an equidistant position from the center 
and the surface of the wire (dashed line). 

the Algi in 1.5 nm SiNWs, restricting ourselves to substitution at the innermost lattice site. 
However, we have found that at such small diameters the dependence is highly non-trivial, 
and a clear trend can hardly be established (see Fig. HI). On the one hand the formation 
energy of the neutral defect varies as fiA\ of the nanoparticle approaches the bulk value; on 
the other hand, the energy at which the Alg ; becomes dominant changes due to the quantum 
confinement effect. 

The impurity concentration Ni is related to the formation energy Ef through Ni = 
N s exp(—Ef/kBT), where N s is the concentration of available sites. We estimate, for T = 
490 °C, an n-type doped Si substrate, and Ef = 0.45 eV as obtained from Fig. @J that the 
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FIG. 4: Formation energy of Al substitutional (Algi) point defect in (110) and (111) SiNWs with 
a diameter of 1.0 and 1.5 nm. The Alsi is at the center of the wire (continuous line in Fig. [3]). The 
electron chemical potential [i e spans the NW band-gap, which is larger for the thinner wires. 

Al concentration for a 1.5 nm (111) SiNW is 5.33xl0 19 cm -3 . This is in agreement with 
the strong Al ^-doping of the SiNWs reported by Wang et al. 23[]. Unfortunately, they can 
only provide a rather loose upper bound on the Al concentration, 10%, due to the limited 
sensitivity of their experimental setup. We notice that such a concentration alone cannot 
explain the observed tapering of the wires based on a catalyst nanoparticle consumption 
mechanism. Hence, either uncatalyzed deposition of Si on the side of the wires 23j or 



catalyst metal outdiffusion 



29 



30] must be invoked. 



Nevertheless, we expect SiNWs to be able to sustain larger Al concentrations than their 
bulk counterparts. Let us consider the vapor-solid-solid growth of p-type layers on top of 
bulk Si. If the substrate is n-doped, Al is easily incorporated in the beginning of the growth 
process. However, as the impurity concentration increases, a space charge zone will develop, 
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bringing the Fermi level closer to the valence band, and the formation energy of Al point 
defects will become larger. In these conditions we predict an Al concentration of 2.3xl0 17 
cm -3 . In the case of SiNWs this compensation mechanism is not expected to be as efficient 
as in thin film growth, because the Al impurities incorporated in the reduced NW volume will 
not significantly move the Fermi level of the whole system, allowing larger Al incorporation. 

In summary, we have presented a generalized formalism that allows to calculate the for- 
mation energy of a defect in one-dimensional semiconductor systems. We avoid using bulk 
derived quantities and we introduce the unambiguously defined chemical potential of the 
nanowire primitive cell. The Madelung correction is extended to the case of arbitrary cell 
shape and dielectric tensor. We apply this formalism to the study of acceptors in silicon 
nanowires grown catalytically, focusing on the stable Al substitutional defects. The rela- 
tively low formation energy for ^-doping conditions further decreases as the Fermi energy 
moves upwards and it finally becomes negative, thus leading to indiscriminate incorporation 
of Al from the catalyst nanoparticle in absence of compensation mechanisms. We have calcu- 
lated the Al concentration for silicon nanowires grown with the vapor-solid-solid mechanism, 
predicting an Al solubility at least one order of magnitude larger than in bulk. 
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